Detection of variable tunneling rates in silicon quantum dots 
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Reliable detection of single electron tunneling in quantum dots (QD) is paramount to use this category of 
device for quantum information processing. Here, we report charge sensing in a degenerately phosphorus- 
doped silicon QD by means of a capacitively coupled single-electron tunneling device made of the same 
material. Besides accurate counting of tunneling events in the QD, we demonstrate that this architecture can 
be operated to reveal asymmetries in the transport characteristic of the QD. Indeed, the observation of gate 
voltage shifts in the detector's response as the QD bias is changed is an indication of variable tunneling rates. 



Silicon-based quantum dot (QD) architectures^^ ^ for 
quantum computing have recently emerged as an attrac- 
tive possibility in view of their scalability, their compat- 
ibility with the widely accessible complementary metal- 
oxide-scmiconductor technologjSl and their long coher- 
ence time for electron spin states.'^ In these systems, 
readout of qubit states typically requires charge sens- 
ing and spin-to-charge conversion.!^ Nevertheless, charge 
sensing is useful to achieve the few-electron regime in 
the dot which is instrumental to perform spin manip- 
ulations.'' Demonstration of reliable single-electron de- 
tection is, therefore, crucial for the implementation of 
quantum logical operations and future development of 
complex quantum computational schemes. 

Here, we report charge sensing using two capacitively 
coupled single-electron tunneling devices (SET) fabri- 
cated in degenerately phosphorus-doped silicon. One 
SET is used as a QD weakly coupled to two electron 
reservoirs (namely, source and drain) while the remain- 
ing SET is used as a detector of the QD charge state. We 
show that modifications of the drain-source bias voltage 
of the QD can lead to changes in the transparency of its 
tunnel barrier (s) . This results in the observation of a neg- 
ative differential conductance (NDC) in the QD's stabil- 
ity plot as well as in a shift with respect to gate voltage of 
the expected detector's response.l^ These effects may arise 
from some intrinsic characteristics of the material, such 
as defects at the Si/Si02 interface^ and disordered dis- 
tribution of dopants in the nano-structure.^ Ultimately, 
they may limit the ability of detecting electrons and need 
to be taken into account for the reliable design of these 
devices. 

The top inset of Fig. [ija) shows a scanning electron 
micro-graph (SEM) image of a device similar to those 
investigated. The system material is silicon-on-insulator 
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with a concentration of phosphorus atoms in the active 
layer of about 10^^ cm~'^. The two SETs are trench- 
isolated via reactive ion etching so that no current can 
flow between them and coupling is only capacitive. The 
islands of the two SETs are electron-beam-defined by 
patterning constrictions in the nano-structure that act 
as tunnel barriers due to severe carrier depletion caused 
by both trapping at the Si/Si02 interfac^^^ and dopant 
deactivation.!^ It has been previously reported that spe- 
cific fabrication arrangementJi^ should be put in place 
to mitigate the effect of disorder on these devices' char- 
acteristics. Our devices were rather realised to optimize 
coupling effects, hence the observation of asymmetries in 
the transport characteristics is not surprising, as we shall 
discuss next. 

In order to significantly suppress thermal fiuctuations 
and allow Coulomb Blockade (CB) to arise, experiments 
have been carried out in a He'^ cryostat at the base tem- 
perature of 300 mK. The top SET is typically used as 
a QD whose electron number is modified by sweeping 
the voltage of the adjacent gate electrode (Vg)- The av- 
erage current fiowing from source to drain through the 
QD is controlled by acting on the bias voltage iyos)- 
The bottom SET is used as a charge detector whose op- 
erating point is selected to maximize charge sensitivity. 
Since the voltages applied to the QD gate to modify its 
electron number slightly affect the potential of the detec- 
tor, simultaneous sweep of the gate electrode close to the 
detector (C) is performed to compensate for unwanted 
shifts away from optimal bias. The suitable compensa- 
tion ratio between the two voltages is selected according 
to capacitive considerations, as detailed elsewhere.'^^ 

The stability plot for the QD device when the detec- 
tor's electrodes are grounded is shown in Fig. [T]^a). Be- 
sides the expected diamond-shaped characteristics due to 
CB, the presence of a NDC for negative values of Yds 
is clearly seen. This feature's onset voltage is modulated 
by Vg in a similar manner to the CB diamonds edges 
(dotted line in the figure). This implies that the electro- 
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FIG. 1. a) QD stability plot. Differential conductance varies 
from -0.5 fiS to +1.0 fiS in the scale range reported. Negative 
values of the differential conductance define the NDC feature 
at negative Yds- Top inset: SEM image of a representative 
device showing both QD and detector. Scale bar is 200 nm. 
Bottom inset: experimental (solid line) and calculated (dotted 
line) IV characteristic of the QD at Vg=-2.55V. b) Schematic 
energy level diagrams showing the origin of the NDC in the 
QD for increasingly negative Yds- The lettering indicates 
which region of the stability plot in (a) each diagram refers 
to. Solid (dashed) lines in the island represent levels available 
(unavailable) for transport. Transport channels are evenly 
separated by the charging energy of the QD. 



static potential of the dot plays a significant role in the 
mechanism originating the NDC. Furthermore, we do not 
observe any NDC for Vds > ^ which indicates an asym- 
metrical transport mechanism consistent with unequal 
tunneling probabilities through the two tunnel barriers. 
The formation of these barriers is affected by the energy 
spread of interface trap^i^ as well as by the random spa- 
tial distribution of donors in the constricted regions.l^ 
These phenomena are likely to give rise to asymmetri- 
cal bias-dependent tunneling probabilities in our device. 
Nonetheless, variable tunneli ng ra tes in Si nano-wires 
have been previously report ed! ^^ * ^^ * In Fig.[ljb) a possible 



mechanism to explain the origin of the NDC is sketched. 
Firstly, we assume a strong asymmetry in the shape of the 
tunnel barriers. In particular, the source-island tunnel 
barrier (SITB) should show energy-dependent tunneling 
probability whilst the drain-island tunnel barrier (DITB) 
should be nearly energy-independent in the bias range of 
interest. Secondly, it is reasonable to assume that the ca- 
pacitive coupling between the QD and the source-drain 
electrodes affects the dot's energy levels whenever the 
Fermi level in either reservoir is modified via Vug. In 
region a of the stability plot, CB is lifted because the 
Fermi energy in the QD, fiQo, {solid line in Fig. [ijb)) 
falls within the transport window between the Fermi en- 
ergies of source (/Xg) and drain {^d)- For more negative 
Yjjs (region /3) the current increases monotonically be- 
cause an extra level enters the transport window. By 
contrast, in region 7, an increasingly negative Yds pro- 
duces a decrease of current (i.e. a NDC) because the 
tunneling probability is reduced given the concurrent ef- 
fects of the bias-dependent SITB and the QD's poten- 
tial shift (note that we assume the drain electrode to be 
grounded, hence a negative Yds would pull down both fig 
and fJ-QD leaving unchanged). Finally, by increasing 
further | Vds I a third transport channel becomes available 
(region 6); this compensates the increased opacity of the 
SITB and the current recovers the increasing trend. The 
bottom inset of Fig. [ija) compares one of the experimen- 
tal 1-Yds curves with the calculated characteristics for a 
parabolic shape of SITI^^^and an abrupt shape of DITB. 
An analytical approactP^ is used to calculate the current 
through a double tunnel junction system as a function of 
bias voltage. The calculation shows that a NDC arises 
due to the mechanism detailed above. Importantly, the 
value of Yds at the onset of the NDC depends on the 
parabolic coefficients of the SITB as well as on the de- 
gree of coupling between the reservoirs and the QD. In 
the calculation these parameters have been chosen to fit 
the experimental results. Since conduction by thermally 
activated carrier is not taken into account in this model, 
the agreement between experiments and calculation is ex- 
pected to become less good as the bias voltage increases. 
This may explain the discrepancy in the region beyond 
the NDC. 

We now turn to discuss the effect of this tunneling 
probability modulation on the charge sensing mechanism. 
In Fig. 2j the QD and detector currents (solid black and 
squared blue lines, respectively) are simultaneously mon- 
itored as Yg is swept for three values of Yds- In (a), 
V£)s=-2mV which falls within region a of the stability 
plot. It can be noticed how sharp changes on the detec- 
tor's trace correlate with the dips observed in the QD's 
current, as the dashed lines indicate. This is the clear sig- 
nature of the abrupt change in the detector's electrostatic 
potential caused by tunnelling events of single electrons 
in the QD. In (b), V£)s=-6mV and we are addressing 
region /3 of the stability plot. It is of note that now the 
sensing is slightly misaligned with respect to the gate 
voltages at which the current dips. In (c), VDS=-13mV 
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FIG. 2. Simultaneous current response of the QD (smooth 
lines) and the detector (squared lines) for different Yds as the 
QD's gate is swept, a) Vi35=-2mV (region a of the stability 
diagram). Sensing steps occur at the same voltages as the 
CB dips, b) Vj3s=-6mV (region /3 of the stability diagram). 
Sensing steps are shifted with respect to the CB dips, c) 
V£)S=-13mV (region 5 of the stability diagram). Increased 
shift is noticeable. The detector's bias point is chosen to 
maximize charge sensitivity but is not the same for all the 
traces. 



which is in region 5, i.e. beyond the NDC. In this case, 
an even larger sensing shift in gate voltage can be ob- 
served. At this bias, the sensing takes place in closer 
proximity to a current peak than a dip. It is also worth 
noting that the amplitude of the CB oscillations is sig- 
nificantly reduced with respect to region (3, despite the 
increase of |Vds|- In order to improve the signal-to-noise 
ratio, the detector's bias point was not the same for the 
three graphs reported; in particular, it was occasionally 
needed to change the detector's drain-source voltage and 



this explains why Idet does not have always the same 
sign. However, this kind of modification in the detector's 
bias point does not affect its charge sensitivity; therefore 
we can confidently attribute the observed shifts to vari- 
able tunneling rate in the QD, as we shall demonstrate 
next. 

As mentioned earlier, in region a only one transport 
channel at a time falls within the transport window and 
electrons enter the QD with a rate and leave with 
a rate (see Fig. |3][a)). The symmetry of the stabil- 
ity plot for -5mV<VD5' <+5mV allows us to consider 
Td ~ Fs for this bias range. Under this condition, the 
QD's current dips whenever figo is equidistant from the 
two reservoirs' Fermi levels. The average electron num- 
ber, N, in the QD is then given bj^^ 

N = No + n^No + — ^ ^No + l (1) 

I d + I s ^ 

where iVo is the excess number of electrons already in 
the dot and n is the fractional number of an electron 
added to the QD as a result of the opening of the trans- 
port channel(s) between the drain and source electrodes. 
When the two barriers are almost symmetrical (e.g. in 
region a), n is very nearly equal to | and the transition 
from A^o to A'o -l- 1 electrons in the dot is sketched in 
Fig. Jsljb) by the solid line. From the measurements of 
Fig. [2[a) we can conclude that, since the dips of Iqd line 
up with the onset of the sensing steps, these latter occur 
whenever the electron number in the QD is N = Nq + ^. 
However, when the barriers are strongly asymmetrical N 
approaches an integer (see eq. ([T])), precisely Aq -|- 1 (Nq) 
for Frf ^ (Frf <^ Fs); as a consequence, the sensing 
points will be shifted towards smaller (larger) gate volt- 
ages. In other words, the detector's response would shift 
in gate voltage on the conditions fj.d = fJ,QD for F^ ^ F^ 
and fig = fJ-QD for ^ Ts- Fig- |3][b) shows with red 
dotted lines the situation relevant to our experimental 
condition. Indeed, the presence of a NDC at Yds < is 
consistent with a strong reduction of Fs for increasingly 
negative bias; this leads to an increase of the fractional 
number of electrons in the QD and the consequent shift 
of the detector's response toward smaller Vq. Further 
confirmation of a strong decrease of tunneling probabil- 
ity comes from the decreasing amplitude of the CB oscil- 
lations for increasingly negative bias voltages (compare 
traces in Fig. [2]jb) and (c)). 

For the sake of completeness, we wish to point out 
that for increasingly negative YdSj apart from making 
the tunnel barriers asymmetrical, more transport chan- 
nels become available and the fractional number of elec- 
trons in the QD will be determined by the tunnelling 
rates relevant to each of these levels. Fig. ^c) repre- 
sents the condition pertinent to region d of the stability 
plot where the deeper energy levels will contribute to the 
higher fractional electron numbers. Indeed, by assuming 
again a parabolic SITB and a sharp DITB, we have 

Firf « F2d « Fsd « Fi, » T2s » F3, (2) 
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FIG. 3. a) QD energy level diagram for the bias region a 
where Td ~ Ts at the center of the transport window, b) the 
effect of sweeping Vg (i.e. the Fermi energy of the QD) over a 
CB oscillation on both the average number of electrons in the 
QD and the current in the detector. Solid lines represent the 
situation Fd ~ Fs (i.e. n;^ |), dotted lines are for Td 3> Fs 
(i.e. n<l). c) QD energy level diagram for the bias region 
S. Different tunneling rate amplitudes are schematically indi- 
cated by the arrow thicknesses. 



^id{s) being the tunneling rate of the i*^ channel to drain 



(source). This would lead to 
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< n2 < na < 1 



(3) 



Hi being the fractional electron number of the i*^ chan- 
nel. Clearly, the overall effect in terms of detector's re- 
sponse shift would be comparable to the one described 
for the single transport channel. 

In conclusion, we have reported charge sensing of a QD 
by a capacitively coupled single-electron transistor. The 
QD has revealed an asymmetrical characteristic such as 
a NDC for negative bias. We have also observed a shift 
in gate voltage of the detector's response while varying 
the QD source-drain voltage from the blockaded region 
towards the NDC region. We have illustrated that both 
effects are consistent with variable tunnelling rates in one 



of the QD's tunnel barriers. The ability to detect mod- 
ifications in the tunnelling probability of a QD shown 
here is to be mastered in order to correctly operate these 
devices for the readout stage in quantum computational 
schemes. 
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